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Abstract 

The  desired  performance  of  a  proton  exchange  membrane  fuel  cell  (PEMFC )  stack  is  always  hindered  by  improper  matching  of  critical  design  and 
operating  parameters  associated  with  the  stack.  This  paper  presents  typical  design  considerations  for  a  small  stack,  which  essentially  operates  under 
non-pressurized  and  non-humidified  conditions.  With  the  objective  of  maximized  energy  efficiency  (or  fuel  economy),  the  relationships  between 
conductivity  and  water  loading,  water  loading  and  relative  humidity,  and  relative  humidity  and  air  stoichiometric  number  were  established  under 
different  temperatures.  Through  these  relationships,  the  conductivity,  which  is  performance  related,  is  correlated  with  the  air  stoichiometric  number. 
Thus,  the  latter  has  become  the  sole  control  parameter  for  the  small  PEMFC  stack  at  a  fixed  temperature  operating  under  the  above  mentioned 
constraints.  The  results  show  that,  under  a  low  air  stoichiometric  number,  the  higher  the  operating  temperature  is,  the  higher  the  conductivity  will 
be.  However,  opposite  trend  can  be  seen  for  high  air  stoichiometric  number.  For  a  fixed  operating  temperature,  the  conductivity  drops  dramatically 
when  the  air  stoichiometric  number  exceeds  the  threshold  value.  The  higher  the  operating  temperature,  the  more  rapid  drop  of  the  conductivity 
against  the  increased  air  stoichiometric  number  will  be. 

©  2005  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Fuel  cell  systems  produce  electrical  energy  at  very  high 
efficiency.  Generally,  the  efficiency  is  better  than  that  of  gen¬ 
erators  driven  by  internal  combustion  engines  such  as  piston 
engines  and  gas  turbines.  Only  large  combined  cycle  power 
plants  can  outperform  fuel  cells.  However,  if  high  electrical 
efficiency  is  required,  fuel  cells  can  be  combined  with  gas 
turbines,  thereby  surpassing  the  efficiency  of  combined  cycle 
systems,  while  maintaining  zero/low  emission  characteristics. 
In  particular,  due  to  the  synergistic  effects  of  the  integrated 
solid  oxide  fuel  cell  and  gas  turbine  technologies,  simulation 
results  showed  that  an  overall  system  efficiency  of  70%  (net 
AC  power  to  LHV  of  fuel,  typically  natural  gas)  or  higher  is 
possible  with  a  more  complex  thermodynamic  cycle  [  1  ] .  In  addi¬ 
tion,  due  to  the  ability  to  integrate  power  production  in  dwelling 
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areas,  efficient  use  of  the  waste  heat  for  various  applications 
is  possible. 

It  is  generally  agreed  that  proton  exchange  membrane  fuel 
cell  (PEMFC)  is  the  most  promising  fuel  cell  technology  for 
mobile  and  portable  applications  due  to  its  high  energy  den¬ 
sity,  low  weight  and  volume  compared  with  other  types  of 
fuel  cells.  In  addition,  its  operating  temperature  of  less  than 
80  °C  makes  it  suitable  for  rapid  start  up.  Literature  review  on 
the  design  of  PEMFC  showed  that  researchers  either  focus  on 
detailed  component  and/or  small  stack  design  with  the  help  of 
sophisticated  models  or  by  tedious  experimental  trial-and-error 
involving  different  conceptual  design.  Lee  and  Lalk  [2]  claimed 
to  have  developed  a  technique,  which  was  specially  designed 
so  that  models  developed  based  on  this  technique  can  be  used 
to  determine  the  fundamental  thermal-physical  behavior  of  a 
fuel  cell  stack  under  any  operating  and  design  parameters.  It  can 
also  be  used  to  replicate  the  performance  of  a  specific  fuel  cell 
system.  Mann  et  al.  [3]  developed  a  generic  model  that  could 
accept  as  input  not  only  values  of  the  operating  variables  such 
as  anode  and  cathode  feed  gas,  pressure  and  compositions,  cell 
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temperature  and  current  density,  but  also  cell  parameters  includ¬ 
ing  active  area  and  membrane  thickness.  A  special  feature  of  this 
model  is  the  addition  of  a  term  to  account  for  membrane  aging 
by  considering  the  deterioration  of  water-carrying  capacity  with 
time.  Jiang  and  Chu  [4]  reviewed  and  evaluated  the  design  of 
three  types  of  PEMFC  structure  under  different  humidity  and 
temperature,  including  bipolar,  pseudo  bipolar  and  monopo¬ 
lar  stacks  and  concluded  that  they  are  suitable  for  moderate 
to  high  power,  high  power  and  low  power  (in  particular  with 
high  voltage)  applications,  respectively.  While  sophisticated  fuel 
cell  models  could  give  a  detailed  description  of  the  fuel  cell 
characteristics  and  performance  subject  to  the  change  of  per¬ 
formance  related  parameters  or  operating  conditions,  simple 
catalyst  model  developed  by  Chan  and  Tun  [5]  could  provide 
the  information  on  the  relative  effect  of  each  key  parameters 
on  the  fuel  cell  performance  so  that  effort  for  improvement  can 
be  focused  on.  Rajalakshmi  et  al.  [6]  developed  a  technique 
for  determining  the  distribution  of  current  density  in  an  operat¬ 
ing  PEMFC  using  segmented  cell  structure.  Real-time  current 
density  distribution  was  measured  which  can  contribute  to  an  in- 
depth  understanding  of  the  reactant  distribution  over  the  active 
fuel  cell  area  for  optimization  of  the  fuel  cell  performance. 
Kumar  and  Reddy  [7]  focused  on  the  improvement  of  PEMFC 
performance  through  optimization  of  the  channel  dimensions 
and  shape  in  the  flow-field  bipolar/end  plates.  The  authors  found 
that  for  hydrogen  consumption  greater  than  80%,  the  optimum 
dimension  values  for  channel  width,  land  width  and  channel 
depth  are  close  to  1.5,  0.5  and  1.5  mm,  respectively.  Their  study 
also  showed  that  triangular  and  hemispherical  shaped  cross- 
section  channel  performed  better  than  square  cross-section  chan¬ 
nel  in  terms  of  hydrogen  consumption.  Chan  et  al.  [8]  developed 
a  PEMFC  model  with  anode  CO  kinetics,  which  is  essentially 
a  model  that  marrying  the  work  of  Bernardi  and  Verbrugge 
[9]  with  that  of  Springer  et  al.  [10].  The  model  estimates  the 
superficial  water  velocity  in  the  membrane  electrode  assembly 
(MEA)  as  well  as  the  poisoning  effect  of  CO  concentration  on 
the  PEMFC  performance.  Hottinen  et  al.  [11]  studied  a  planar 
cell  design  for  free-breathing  fuel  cell,  in  which  the  cathode 
side  of  the  cell  is  directly  exposed  to  the  air,  while  the  anode 
side  flow-field  geometry  is  of  ‘Z’  shape.  The  authors  claimed 
that  this  kind  of  design  is  feasible  for  small-scale  applications, 
such  as  portable  computers.  Carnes  and  Djilali  [12]  presented 
an  algorithm  for  estimation  of  PEMFC  model  parameters  using 
a  constrained  nonlinear  least  square  algorithm.  Estimation  of 
five  model  parameters  (i.e.,  membrane  conductivity,  anode  and 
cathode  current  density,  oxygen  diffusion  coefficient  in  gas  dif¬ 
fusion  layer  and  catalyst  layer)  in  a  simple  ID  electrochemistry 
model  using  two  different  experimental  polarization  curves  has 
been  demonstrated  using  this  procedure. 

Whatever  approach  the  researchers  adopted,  the  ultimate  goal 
is  to  achieve  improved  PEMFC  performance  through  optimiza¬ 
tion  of  key  design  parameters  under  the  constraint  of  opera¬ 
tional  requirements.  It  seems  that  between  tedious  experimen¬ 
tal  trial-and-error  and  sophisticated  modelling,  there  is  a  need 
to  establish  relationships  for  matching  of  critical  parameters 
in  the  PEMFC  stack.  In  this  paper,  the  authors  present  typi¬ 
cal  design  considerations  for  a  small  stack,  which  essentially 


operates  under  non-pressurized  and  non-humidified  conditions. 
With  the  ultimate  goal  of  maximized  energy  efficiency  (or  fuel 
economy),  the  relationships  between  conductivity  and  water 
loading,  water  loading  and  relative  humidity,  and  relative  humid¬ 
ity  and  air  stoichiometric  number  under  different  constant  tem¬ 
peratures  have  been  established.  Through  these  relationships,  the 
conductivity,  which  is  performance  related,  is  correlated  with  the 
air  stoichiometric  number.  Thus,  the  latter  has  become  the  sole 
control  parameter  for  the  small  PEMFC  stack  at  a  fixed  tem¬ 
perature  operating  under  the  above  mentioned  constraints.  Note 
that,  since  the  focus  is  placed  on  small  PEMFC  stacks,  the  tem¬ 
perature  range  of  25-40  °C  under  study  should  be  reasonable. 

2.  Non-pressurized,  non-humidified,  dead-ended  anode 
flow  channel  PEMFC 

A  fuel  cell  that  is  an  ambient  pressure  air-breather  is  designed 
to  use  low-pressure  hydrogen  on  the  anode  side  of  the  membrane 
electrode  assembly  and  ambient  air  on  the  cathode  side.  This 
type  of  the  cell  can  be  designed  to  have  a  dead-ended  anode 
flow  channel  with  a  pre-set  purge  cycle,  which  is  sometimes 
referred  to  as  burping  the  cell.  A  valve  may  be  used  to  release  a 
small  amount  of  hydrogen  from  the  anode  flow-field  in  the  other 
side  of  the  flow  channel  where  condensed  liquid  water  and  trace 
impurities  of  the  fuel  trapped  in  the  anode  flow  channel  could 
be  scavenged.  In  this  case,  it  can  be  assumed  that  all  hydrogen 
introduced  to  the  cell  will  be  consumed  and  this  consumption 
can  be  determined  from  the  output  current  of  the  cell. 

Water  plays  a  critical  role  in  proton  exchange  membrane  fuel 
cells.  Both  the  Nation  membrane  and  the  Nation  material  in  the 
catalyst  layers  need  to  be  hydrated  in  order  to  conduct  protons 
efficiently.  The  protonic  conductivity  increases  with  the  increase 
of  water  content.  To  achieve  enough  hydration,  water  is  normally 
introduced  into  the  cell  externally  by  a  variety  of  methods  such  as 
liquid  injection,  steam  introduction  and  humidification  of  reac¬ 
tants  by  passing  the  hydrogen  and  air  through  humidifiers  before 
entering  the  cell.  Humidification  by  the  last  method  is  relatively 
easier  to  handle  and  therefore,  it  is  the  most  commonly  used 
technique.  To  simplify  the  system  and  to  limit  the  amount  of 
water  to  be  carried  with  the  fuel  cell  for  the  humidification  pur¬ 
pose,  the  water  generated  by  the  fuel  cell  at  the  cathode  should 
be  recycled  as  much  as  possible. 

The  design  of  a  fuel  supply  system  for  non-humidified  and 
non-pressurised  fuel  cell  is  less  costly.  Since  the  cell  is  dead- 
ended,  there  is  no  need  to  meter  or  regulate  the  fuel  flow  rate. 
This  means  that  no  mass  flow  controllers  are  required,  except 
one  that  is  needed  for  measuring  the  fuel  consumption.  The 
cell  consumes  fuel  as  required  from  a  source  of  low-pressure 
hydrogen.  If  there  is  no  load  connected  on  the  cell,  the  fuel 
consumption  is  close  to  zero.  This  type  of  cell  is  self-heating  up 
to  the  equilibrium  temperature  and  is  suitable  for  power  output 
smaller  than  100  W. 

3.  Model  development 

The  protonic  conductivity  of  a  polymeric  membrane  is 
strongly  dependent  on  membrane  structure  and  membrane  water 
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Water  production 


Fig.  1.  Schematic  of  water  sources  and  water  fluxes. 


loading  or  membrane  hydration  (1).  X  is  defined  as  the  number 
of  water  molecules  per  sulfonate  acid  groups. 

,  =  iV(H20) 

/V(S03H) 

Thus,  a  key  challenge  in  the  evaluation  of  ionomeric  mem¬ 
branes  for  fuel  cell  has  been  focused  on  the  analysis  of  combined 
structural  and  water  uptake  characteristics  required  to  achieve 
the  highest  protonic  conductivity. 

Sources  of  water  and  the  various  flux  components,  which 
redistribute  water  in  an  operating  PEMFC,  are  shown  schemat¬ 
ically  in  Fig.  1 .  Water  is  produced  at  the  cathode  by  the  elec¬ 
trochemical  process  and  enters  the  cathode  via  the  wet  air  as 
water  carrier.  Water  is  transported  through  the  membrane  from 
the  anode  to  the  cathode  by  electro-osmotic  drag,  i.e.,  by  a  water 
flux  associated  with  the  protonic  current.  The  flux  of  water  due 
to  electro-osmotic  drag  (mol  cm-2  s-1)  is  given  as: 

_  -fr?drag(k) 

^vw,drag  —  p 

where  I  is  the  cell  current  (A  or  C  s-1);  n^rag  (A.),  the  electro- 
osmotic  drag  coefficient  (cm-2)  at  a  given  state  of  membrane 
hydration  X\  and  Tis  the  Faraday  constant  (C  mol-1).  This  water 
migratory  flux  tends  to  dry  out  the  anode  side  of  a  cell  under 
current  and  adds  to  the  faradaic  production  of  water  at  the  cath¬ 
ode,  thus  producing  excess  water  at  the  cathode.  This  build-up 
of  water  at  the  cathode  is  relieved  by  diffusion  down  the  result¬ 
ing  water  concentration  gradient  and  by  hydraulic  permeation  of 
water  in  differentially  pressurised  cells  where  the  cathode  is  held 
at  higher  overall  pressure.  The  fluxes  (mol  cm-2  s-1)  brought 
about  by  the  latter  two  processes  within  the  membrane  are: 

(Vw.diff  =  —D(X)  — 

A  z 


N w.hyd  —  &hyd(A) 


A p 
A  z 


where  D  is  the  diffusion  coefficient  (cm2s-1)  in  ionomer 
of  water  content  A;  Ac/ A",  the  water  concentration  gradient 
(mol  cm-4)  along  the  "-direction  of  membrane  thickness;  &hyd> 
the  hydraulic  permeability  (molsg-1)  of  the  membrane;  and 
A  pi  A"  is  the  pressure  gradient  (g  cm-2  s-2)  along  z.  The  steady- 
state  water  profile  across  the  ionomeric  membrane  for  given  cell 
current  density,  external  humidification  conditions  and  differen¬ 
tial  pressurisation,  is  the  resultant  of  electro-osmotic,  diffusive 
and  hydraulic  fluxes. 

Modelling  work  on  the  PEM  fuel  cell  began  as  electrode- 
membrane-electrode  models.  Most  previous  models  focused  on 
the  fully  hydrated  membrane  fuel  cell  with  the  condition  that 
reactants  are  thoroughly  humidified  prior  to  entering  into  the 
fuel  cell.  The  membrane  is  liquid  water-equilibrated  membrane 
and  the  electrodes  are  hydrophobic  in  nature,  which  are  made  of 
the  catalyst  mixed  with  PTFE,  then  brushing  liquid  membrane 
solution  onto  the  catalysed  face  of  the  electrode  and  allowing  the 
liquid  to  seep  into  the  electrode’s  pores.  In  this  case,  the  conduc¬ 
tivity  is  only  a  function  of  temperature.  For  example,  the  conduc¬ 
tivity  ofNafionis  0.17  S  cm-1  at80°Cand0.21  Scm-1  at95°C 
[9].  However,  for  fuel  cells  that  operate  at  low  temperatures 
and  reactants  are  not  humidified  or  only  slightly  humidified, 
the  membrane  is  water  vapour-equilibrated  membrane  and  the 
electrodes  are  usually  hydrophilic  ones,  which  are  made  of  the 
catalyst  mixed  with  liquid  membrane  solution.  The  behaviour 
of  water  vapour-equilibrated  membrane  is  much  more  complex 
than  liquid  water-equilibrated  membrane.  In  this  case,  the  con¬ 
ductivity  is  a  function  of  the  operating  temperature  and  airflow 
rate.  The  latter  is  applied  to  this  study. 

The  conductivity  of  membrane  is  a  function  of  water  loading. 
If  water  loading  is  less  than  two,  the  Nation  membrane  behaves 
like  an  insulator.  Above  this  threshold,  the  conductivity  rises 
dramatically  with  water  loading,  reaching  a  plateau  in  the  semi¬ 
conductor  range  of  the  order  of  about  10-1  Scm-1  [13].  The 
water  loading  is,  in  turn,  related  to  the  activity  of  water  (or  the 
relative  humidity),  while  the  latter  is,  in  turn,  determined  by  the 
amount  of  water  produced  in  reactions  and  the  humidity  of  the 
air  carrying  to  cathode.  The  amount  of  the  air  carrying  to  the 
cathode,  in  this  study,  is  expressed  as  a  function  of  air  stoichio¬ 
metric  number  and  is  defined  as  molar  ratio  of  the  actual  inlet 
air  supply  to  stoichiometric  air  in  the  hF-air  reaction. 

Since  anode  is  dead-ended  in  this  study,  the  air  stoichiometric 
number  and  fuel  cell  operating  temperature  are  the  only  two 
controllable  parameters. 

Fig.  2  shows  the  interdependency  between  any  two  key 
parameters  that  are  related  to  the  performance  of  PEMFC.  Once 
each  dependency  is  established,  the  relationship  between  the 
conductivity  and  the  air  stoichiometric  number  under  a  fixed 
operating  temperature  can  be  obtained. 


Fig.  2.  Interdependency  between  key  parameters  under  a  fixed  temperature. 
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4.  Relationship  between  conductivity  and  water  loading 


For  the  case  of  a  proton  exchange  membrane  consisting  of 
water  as  the  solvent  and  hydronium  ion  as  the  protonated  solvent, 
the  conductivity  model  developed  by  the  Thampan  et  al.  [13]  is 
expressed  as: 


<r=(e-  s0)9 


cha,o« 


where  s  is  volume  fraction  of  water  in  hydrated  membrane; 
or  wet  porosity,  eo,  the  percolation  threshold  volume  fraction 
of  water  in  hydrated  membrane,  q,  the  Bruggemann  or  critical 
exponent;  X® ,  equivalent  conductance  for  H3O-  at  infinite  dilu¬ 
tion;  (5,  ratio  of  mutual  to  matrix  effective  diffusion  coefficients; 
c'ha.o,  concentration  of  membrane  acid  groups;  and  a  is  the 
degree  of  acid-group  dissociation.  And 

X 

£=  Vu/Vi  +  k 


CHA.O  = 


A .V2 


where  X  is  water  loading,  or  the  number  of  water  molecules  per 
— SO3H  group;  VM,  the  partial  molar  volume  of  membrane;  and 
V2  is  the  partial  molar  volume  of  water,  V?  =  18  cm3  mol-1 . 

EW 

Vm  ~  - 

Po 

where  EW  is  membrane  equivalent  weight,  for  Nafion  117, 
115,  1135  and  112,  it  is  1100;  po  is  density  of  dry  membrane, 
po  =  2.05  gem-3. 


—  ^1,298  eXP 


where  A®  29&  is  equivalent  conductance  for  H3O-  at  infinite 
dilution  at  298  K,  A®  29g  =  349.8  S  cm-2  eq_1 .  E ^  is  activation 
energy  for  viscosity  of  water,  =  14  kJ  mol“ 1 .  R  is  universal 
gas  constant,  R  =  8.3143  J  mol-1  K.  T  is  the  operating  tempera¬ 
ture. 


(A  +  1)  -  \J (X+  l)2  —  4A(1  —  I/A'a.c) 
2(1  -  1/^a.c) 


where  K\  q  is  equilibrium  constant  for  proton  solvation  in  terms 
of  concentrations. 


Aa.c  =  Af\,c,298exP 


A  H° 

(l 

1 

R 

\T 

298  ) 

where  A'a.c, 298  is  equilibrium  constant  for  proton  solvation  in 
terms  of  concentrations  at  298  K,  A'a.c.298  =  6.2.  A H°  is  the 
enthalpy  change,  A H°  =  —52.3  kJ  mol-1 . 

The  critical  exponent  q  is  universal  constant  predicted  to  be 
about  1 .5.  The  threshold  value  eo  is  best  determined  from  exper¬ 
iments  as  a  fitting  parameter  with  value  of  0.06.  <5 ,  a  value  of  5.5, 
is  judiciously  adjusted  according  to  the  experiment. 

Fig.  3  shows  the  relationship  between  the  conductivity  of  the 
membrane  and  the  water  loading  (the  number  of  water  molecules 


Fig.  3.  Relationship  between  conductivity  of  the  membrane  and  water  loading. 


per  — SO3H  group)  under  different  operating  temperature.  It  is 
noted  that  when  the  water  loading  is  below  five,  the  conductivity 
drops  significantly.  For  a  fixed  water  loading,  the  conductivity 
improves  with  the  increase  of  temperature.  The  effect  of  temper¬ 
ature  on  conductivity  diminished  when  the  water  loading  drops 
toward  the  threshold  value  of  two. 

5.  Relationship  between  water  loading  and  relative 
humidity 

Membrane  water  content/loading  (A)  is  essentially  the  ratio 
of  the  number  of  water  molecules  per  SC>3_H+  charge  site.  It 
can  be  determined  from  the  activity  of  water  vapour  at  the  elec¬ 
trode/membrane  interface.  The  activity  (a 2)  in  the  vapour  phase 
is  defined  as: 


where  (p  is  the  relative  humidity,  pw  and  /jsat  are  the  pressure 
and  saturation  pressure  of  the  water  vapour,  respectively.  The 
latter  is  fitted  to  the  following  empirical  expression  using  the 
thermodynamic  data  of  the  water: 

log10psat  =  —2.1794  +  0.02953T  —  9.1837  x  10“5T2 
+  1.4454  x  10“7r3 


Zawodzinski  et  al.  [14]  measured  and  expressed  the  water 
loading  as  a  function  of  water  activity  for  the  Nafion  117  mem¬ 
brane.  Fig.  4  in  ref.  [14]  shows  the  experimental  fit  of  “isopiestic 
curve”  measured  at  30  °C.  This  experimental  fit  of  A  versus  a 
can  be  expressed  by  the  following  equation. 

A(30C)  =  0.043  +  17.8 1«2  —  39.85a2  +  36a3  for  0  <  a2  <  1 

When  a2  =  1 ,  the  value  of  A  in  equilibrium  with  saturated  water 
vapour  is  14  (or  more  specifically,  14.003)  FEO  molecules  per 
charged  site.  When  the  membrane,  following  partial  drying,  is 
immersed  in  liquid  water  at  80  °C,  A  can  go  as  high  as  16.8.  The 
A  is  assumed  to  increase  linearly  from  14  to  16.8  when  the  mole 
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fraction  of  water  exceeds  the  saturation  from  xw,sat  to  3xw.sat, 
i.e.. 


X  =  14  +  1.4(02  —  1)  for  1  <  «2  <  3 


Alternatively,  the  water  sorption  characteristics  of  Nafion  can 
be  well  modelled  by  Brunauer-Emmett-Teller  (BET)  equation 

[15]: 

X  [Ca2/(1  -  o2)][l  -  (H2  +  l)fl^2+1  +  n2af+l] 

Xm  “  1  +  (C  -  l)o2  -  Ca"2+1 

where  A.m  is  the  water  loading  at  monolayer  coverage,  C  is  BET 
constant,  02  is  the  activity,  or  relative  humidity,  and  n2  is  total 
number  of  water  layers  sorbed  on  the  pore  surface. 

The  monolayer  coverage  /,m  was  estimated  from  knowledge 
of  the  specific  pore  surface  area  S  by: 

5  EW 

3-m  —  77  r~ 

Po  N\A2 


where  the  surface  area  occupied  by  an  adsorbate  molecule  on 
the  pore  surface  and  can  be  estimated  from: 


A2 


1.091 


mw2  \  2/3 

PiNk) 


For  Nafion,  S  is  210  m2  cm-3.  /V,\  is  the  Avogadro’s  constant; 
MW2,  the  molecular  mass  of  water;  and  pi  is  the  density  of  water. 
C  and  «2  are  treated  as  fitting  parameters,  C  =  150  and  n2  =  13.5. 

Fig.  4  shows  the  adsorption  isotherm  for  water  loading  by 
Nafion  equilibrated  with  water  vapour.  If  water  loading  of  five  is 
required,  the  relative  humidity  should  be  over  60%.  Combining 
the  relationship  between  conductivity  and  water  loading  and  the 
relationship  between  water  loading  and  relative  humidity,  the 
relationship  between  conductivity  and  relative  humidity  can  be 
obtained,  which  is  shown  in  Fig.  5. 


Fig.  4.  Adsorption  isotherm  for  water  loading  by  Nafion  from  water  vapor. 


Fig.  5.  Relationship  between  conductivity  of  membrane  and  relative  humidity. 


6.  Relationship  between  relative  humidity  and 
stoichiometric  number  of  air 

The  water  activity  is  contributed  from  both  the  wet  air  and 
water  arising  from  the  electrochemical  reaction  of  H2— O2.  If 
the  entire  system  is  taken  to  be  held  at  a  constant  temperature, 
and  the  water  vapour  is  assumed  to  be  an  ideal  gas  and  well 
mixed  in  the  cathodic  chamber,  then  the  partial  pressure  of  the 
water  can  be  calculated  from  the  known  conditions  of  the  air 
(relative  humidity  and  temperature  of  the  atmospheric  air)  and 
the  operating  air  stoichiometric  number. 

The  overall  reaction  can  be  expressed  as: 

mH2  +  v(Oi  +  /6N2)  +  ixTEO  =  aTEO  +  \02  +  zNE 

StQ&npjf 
2(Pt  -  </>in/4at) 

_  StQ  fan  Pig  _  StO  p,  1 

2(pt  -  <PinP\nY  y  2(pt  —  4>mP^)  2’ 

/IS  to  pt 

2(pt  -  (pinPin) 

where  Sto  is  the  stoichiometric  number  of  air;  pt  the  total  pres¬ 
sure;  (p\n  the  relative  humidity  of  inlet  air,  which  is  assumed  to  be 
70%  in  this  study;  the  saturated  vapour  pressure  evaluated 
at  25  °C;  and  /l  =  3.762,  which  is  the  volume  ratio  of  nitrogen  to 
oxygen  in  the  air. 

Thus,  water  partial  pressure  is: 

Pw  _  2(pt  ~  (pinPln)  +  StO</>in 

Pt  Pt  +  (1  +  y6)StO  Pt  +  (pin Pin  ( StO  -  1) 

Fig.  6  shows  the  relationship  between  stoichiometric  number 
of  the  air  and  relative  humidity.  With  the  same  stoichiometric 
number,  the  higher  the  operating  temperature,  the  lower  the  rel¬ 
ative  humidity  will  be.  If  the  stoichiometric  number  is  less  than 


In  the  equation, 


u  =  1 ,  i>  = 


Sto  pt 


2(pt-<pmP^y 


and 
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Fig.  6.  Relationship  between  stoichiometric  number  of  the  air  and  relative 
humidity. 


T  (K) 

Fig.  8.  Relationship  between  conductivity  of  membrane  and  operating  temper¬ 
ature. 


a  critical  value  (i.e.,  14.7,  10.5,  7.7,  and  5.8  for  temperature  of 
25,  30,  35  and  40  °C,  respectively),  the  relative  humidity  can 
always  reach  100%. 

7.  Relationship  between  conductivity  and  stoichiometric 
number  of  air 

Finally,  the  relationship  between  conductivity  of  the  mem¬ 
brane  and  stoichiometric  number  of  the  air  is  established,  which 
is  shown  in  Figs.  7  and  8.  In  Fig.  7,  one  can  see  that  when  the 
air  stoichiometric  number  is  low,  the  conductivity  will  remain 
the  same  over  a  range  of  air  stoichiometric  number  because 
the  membrane  is  wet  under  the  relative  humidity  of  100%.  The 
results  show  that  when  the  stoichiometric  number  is  increased 
above  a  certain  threshold,  the  conductivity  decreases.  It  is  also 
shown  that  the  higher  the  operating  temperature,  the  higher 
the  decreasing  rate  of  conductivity  will  be.  In  Fig.  8,  in  the 
temperature  range  of  298-313  K,  one  can  see  that  at  low  air  sto¬ 
ichiometric  number,  the  higher  the  operating  temperature,  the 


Fig.  7.  Relationship  between  conductivity  of  membrane  and  stoichometric  num¬ 
ber  of  the  air. 


higher  the  conductivity  will  be;  in  the  middle  air  stoichiometric 
number  range,  the  conductivity  increases  with  temperature  until 
reaching  a  peak  then  decreases;  at  high  air  stoichiometric  num¬ 
ber,  the  peak  conductivity  is  vague  and  the  conductivity  mostly 
decreases  with  the  increase  of  operating  temperature. 

8.  Conclusions 

With  the  objective  of  maximized  energy  efficiency,  the  rela¬ 
tionships  between  conductivity  and  water  loading,  water  loading 
and  relative  humidity,  and  relative  humidity  and  air  stoichio¬ 
metric  number  under  different  temperatures  were  established. 
Through  these  relationships,  the  conductivity,  which  is  perfor¬ 
mance  related,  is  correlated  with  the  air  stoichiometric  number. 
Thus,  the  latter  has  become  the  sole  control  parameter  for  the 
small  PEMFC  stack  operates  at  a  fixed  temperature  under  the 
above  mentioned  constraints.  The  results  showed  that,  in  the 
temperature  range  of  298-313  K,  at  low  air  stoichiometric  num¬ 
ber,  the  higher  the  operating  temperature,  the  higher  the  conduc¬ 
tivity  will  be;  in  the  middle  air  stoichiometric  number  range,  the 
conductivity  increases  with  temperature  until  reaching  a  peak 
then  decreases;  at  high  air  stoichiometric  number,  the  peak  con¬ 
ductivity  is  vague  and  the  conductivity  mostly  decreases  with 
the  increase  of  operating  temperature.  For  a  fixed  operating 
temperature,  the  conductivity  drops  dramatically  when  the  air 
stoichiometric  number  exceeds  the  threshold  value.  The  higher 
the  operating  temperature,  the  more  rapid  drop  of  the  conduc¬ 
tivity  against  the  increased  air  stoichiometric  number  will  be. 
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